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We study a quantum computer with fixed and permanent interaction of diagonal type 
| between qubits. It is controlled only by one-qubit quick transformations. It is shown how 

to implement Quantum Fourier Transform and to solve Shroedinger equations with linear 
and quadratic potentials by a quantum computer of such type. The method is adaptable 
| to the wide range of interactions of diagonal form between qubits and to the case when 

C*~) ■ different pairs of qubits interact variously. 



^ ■ 1 Introduction 

o . 

The main difficulty in practical implementation of quantum computing is to fulfill two qubits 
transformations playing a crucial role in quantum algorithms. To perform such transformations 
specially we must in fact artificially and exactly control the degree of their entanglement that 
is determined by overlapping of their spatial wave functions. However to distinguish different 
q-i qubits the share of overlapping amplitudes must be much less than the overall amount of ampli- 

tude and thus in the same degree one-qubits transformations are easier to fulfill than two qubits 
transformations. While we perform a transformation with one pair of qubits a physical interac- 
tion between other pairs of qubits cannot be stopped. This permanent interaction existing in all 
real systems requires special and nontrivial methods of correcting. These difficulties complicate 
a straightforward implementation of quantum algorithms. Here we shall study a nonstandard 
model of quantum computer possessing formally more narrow possibilities but which may be 
more feasible. It is controlled by only one-qubit transformations whereas two qubits interactions 
are fixed and determined by the spatial disposition of qubits, they go permanently in course of 
computation. To show the possibilities of such model we shall at first study how to implement 
Quantum Fourier Transform (QFT) by such quantum computer assuming that the potential of 
two qubits interaction has a diagonal form and decreases as Yukawa potential. Then it will be 
shown how this method can be generalized to wide range of interactions of diagonal form. By 
means of properly chosen one-qubit transformations this method can be easily adopted to the 
case when different pairs of qubits interact differently. At last this approach will be applied to 
the solution of Shoedinger equation for linear and quadratic potentials. 
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Quantum Fourier Transform (QFT) is a key subroutine in quantum computing. It is used in 
variety of algorithms (look at flShl , |AL| , |Ozj| ) as a main step generating interference of amplitudes 
which makes quantum computations so powerful. A simple quantum gate array implementing 
the reversal for QFT is shown at the picture 1. It was proposed in several works and was 
used by Shor for fast factoring (see |{Sh|| ). Let us agree to represent an integer of the form 
a = cio + clq2 + . . . + a\-\l l ~ x by the basic state |ao ai . . . ai-i ) = \a) forming a basis for input 
states of a gate array and dispose all a,j from top to bottom. The same agreement will be for 
output only binary figures bj for an integer b = b + b 2 + . . . + 6;_i2 i_1 will be written in reverse 
order - from bottom to top. 
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Picture 1. Gate array for QFT" 1 with one and two qubits control. Circles 
denote Hadamard gates, two qubits gates has the form (ED) 



This array fulfills QFT" 1 in 0(l 2 ) steps whereas its matrix is N = I 2 dimensioned. 

However, a direct implementation of this scheme requires a control over two qubits trans- 
formations and thus it does not fit into our model of quantum computer. In this paper it is 
shown how QFT and its reversal can be implemented by means of fixed Hamiltonian of two 
particles interaction of the form 
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where all p = po~ — -; b = const; r is a distance between the particles and p\ + p 4 ^ p 2 + P3. 
Dispose our I qubits on one line with equal intervals. We shall consider a case when interaction 
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between jth and kth qubits will have Hamiltonian Hj^ of the forms (0). 

This type of Hamiltonians appears for example in Ising model for particles with spin 1/2. 
The required decrease of interaction with the distance could be obtained by placing each particle 
in the appropriate potential hole. Choosing appropriate unit of the length we can make 6 = 1. 
At first we shall study fixed interaction of the form ([I], A) and then extend our results to ([I], 
B). 



2 Implementation of QFT in within phase shifts 



We assume that QFT and its reversal have the form []: 
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The reversal transformation for QFT can be fulfilled by the following gate array. 
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Picture 2. Frameboxes denote fix interaction of the form (|l], A), 
circles denote Hadamard gates. 



Here by a framebox we denote a unitary transformation of the form U = e tH where 
H — Hj k, and each of Hj & has the form (JT1 , A) with p = ir, r = j — k. If we choose 

l>j>k>0 

a unit of time so that Plank constant multiplied by po equals ir and a unit of length so that 



1 This agreement corresponds to the definition of ordinary Fourier transform. Typically in quantum comput- 
ing literature it is assumed a reversal definition. 
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r = j — k then U will be exactly the transformation of state vector induced by the considered 
Hamiltonian in the unit time frame. We assume here that the time of all one-qubit gates action 
is negligible so that two qubits interaction cannot corrupt phases while these gates act. This 
gate array may be obtained from the previous by insertion of "missing" gates corresponding 
to interactions existing physically in the system with constant Hamiltonian. To prove that 
this gate array fulfils QFT -1 we follow the method of amplitudes counting proposed in the 
paper |JSh| . Given a basic input state \a) consider the corresponding output state. It is a 



linear combination of basic states \b) with some amplitudes. All modules of these amplitudes 
equal 1/y/L and we should only count the phases. For the simplicity we introduce the notation 
a'j = ai_i_j, j = 0, 1, — 1. In the process of the gate array application values of qubits 
with numbers j and k < j pass through the gates on the picture 2 from left to right. Following 
this passage we separate the following four types of segments: interacting of a'j with itself and 
a' k with itself by Hadamard gates, interacting of a'j with a' k (j > k), interacting of a'j with b k for 
j > k, and interacting of bj with bk (j > k). The times of these actions are: zero, k, j — k and 
I — 1 — j correspondingly. Summing the deposits of all these actions we obtain the resulting 
phase 

a j a k k . a!jb k (j-k) bjb k (l - j - 1) 
7i -_ k( ■ T\ L^i -k( ■ T\ ^ ajbj + 7T . , — — . (3) 

i>j>k>o 21 U ~ k ) i>j>k>o 2J U ~ k > i>j>o i>j>k>o 2J U - k ) 

Denote the first and last summands by A and B correspondingly. Their deposits correspond 
to the actions of diagonal Hamiltonians on \a) and \b) correspondingly. Leave this deposit so 
far - till the next section. Take up a part of sum formed by the second and third summands. 
After the replacement j by I — 1 — j this part acquires the form 

tt E ^^ + vr E a*-i-A =2^ E a -^ = 2nS + 27r £ ^11 = 

l-l>k+j>0 2-l>J>0 l>k+j>0 l>k,j>0 

27TS + 27lf 

(4) 

for some integer S. The first summand here does not change the phase and we obtain all what 
is required for QFT -1 but deposits of A and B. 



3 Correcting of phase shifts 

To cope with the deposit of diagonal summands A and B to the phase we present one trick. At 
first consider only a summand A. It consists of addends of the form A,- fc = Cj^a'jO! k1 where Cj tk 
depends only on j and k but not on a. Declare jth and kth qubits separated. We shall apply 
one qubit gate NOT several times to all qubits but separated ones to suppress all two qubits 
interactions excluding interaction between separated qubits. 

At first consider a pair of not separated qubits with numbers p,q, q > p. Their permanent 
interaction in time frame At gives the addend d Piq At a' p a' q to the phase where a real number 
d p>q depends only on how fast does interaction decrease and not on a' p , a' . For example for 
a decreasing of Yukawa type we have d p>q = e~' 9 ~ p '/|g — p\. Now invert one of these qubits, 
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no matter which, say qth by NOT gate. Its state will be 1 — a'. Then the second At period 
of permanent interaction gives the addend d p<q At a' p (l — a! q ) to the phase. At last restore the 
contents of qth qubit by the second application of NOT. The resulting phase shift of this four 
steps transformation will be d Ptq At a' p and it depends on the contents of pth qubit only. Now 
we can compensate this phase shift by a simple one-qubit transformation. If we consider a 
pair of qubits with numbers p, q where one, say pth is separated and other is not, then we can 
compensate their interaction by the same way using one-qubits operations: two NOTs for q 
and one phase shift for pth. 

Now we should so modify this method that compensate all influence of not separated qubits 
simultaneously. For each not separated qubit number p consider the Poisson random process 
A p generating time instants < t \ < t\ < ■ ■ ■ < t v mj> < 1 with some fixed density A > 1, 
Let all A p are independent. Now fulfil transformations NOT on each pth qubit in instants 
t p m sequentially. In instant 1 fulfil NOT on pth qubit if and only if m p is odd. Thus after 
this procedure each qubit restores its initial value. Count the phase shift generated by this 
procedure. Interaction between separated qubits will be unchanged. Fix some not separated 
qubit number p and count its deposit to phase. It consists of two summands: the first comes 
from interaction with separated and the second - from interaction with not separated qubits. 
Count them sequentially. 

1. In view of big density A of Poisson process A p about half of time our pth qubit will be 
in state a' p and the rest half - in 1 — a' p . Its interaction with a separated qubit, say jth brings 
the deposit \d P) ja! p d^ + \d P j{\ — a' p )a'j that is \d p ^a'y 

2. Consider a different not separated qubit number q ^ p. In view of independence of time 
instants when NOTs are fulfilled on pth and qth. qubits and big density A these qubits will be 
in each of states (a' p , a'), (a' p , 1 — a' q ), (1 — a' p , a'), (1 — a' p , 1 — a') approximately a quarter 
of time. The resulting deposit will be \d p , q [a' p a' q + a' p (l — a') + (1 — a' p )a' q + (1 — a' p )(l — a' )] 

4 u p,q- 

A total phase shift issued from the presence of not separated qubits in our procedure now 
is obtained by summing values from items 1 and 2 for all p ^ {j, k}. It is 

\\ 12 d Pd a 'j+ 12 d p^k] + \ 12 d p,i- 
p${j,k} P<£{jM p,q£{j,k} 

This shift can be compensated by one-qubit operations because the first two summands depend 
linearly on the qubits values and the second does not depend on qubits values at all. Thus we 
obtain a gate with permanent two qubits interaction and one-qubit operations fulfilling phase 
shift to dj^djd'k- If we take time frame At instead of unit time in this procedure we obtain 
the phase shift to At dj^d^d k . If we want to obtain the shift to —At dj^a'jd'k we should at 
first apply NOT to the jth qubit, then apply the above procedure, then again apply NOT to 
jth qubit and at last add —At dj^a' k by one-qubit operation. So we are able to make any 
addition of the form c • a'ja' k to the phase for real c independently of its sign. An appropriate 
combinations of such gates gives us a phase shifts to 

12 c i^jd k (5) 
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for any Cj^. Placing such gates before and after QFT - procedure from the previous section 
we compensate summands A and B in the phase and thus obtain a gate implementing QFT -1 . 

Errors arising in this method issue from the possible imperfection of Poisson processes 
generating instants for the qubits inversions and interactions continuing in course of these 
inversions. It can be minimized by increasing of density A and by decreasing of the time needed 
for NOT operation in comparison with typical time of two qubits operation defined by the 
interaction dj t k- 

4 Possibility to use interactions of the different types 

Up to now we considered Hamiltonian of two qubits interaction given by (JTJ, A) with decreasing 
of Yukawa type. Now we are going to extend our technique at first to Hamiltonians of the form 
(H|, A) for arbitrary degree of decreasing, then to Hamiltonians of the form ([[], B) for arbitrary 
decreasing (even for different degrees of decreasing for different pairs of qubits), at last - to the 
interactions which can be diagonalized by one-qubit transformations. Namely, we shall prove 
that it is possible to generate a phase shift of the form (|^) and of the form 

c j,ka'jh (6) 

l>j>k>0 

for any c jjk - 

Consider an interaction of the form (fj, A) for arbitrary degree of decreasing. Generating of 
quadratic phase shifts of the form ([5|) remains unchanged. The only thing we should do is to 
generate QFT -1 in within quadratic phase shift. To do it we start with a gate array represented 
at the picture 2 but now a time interval AT between nearest Hadamard gates will not equal 
1 - it will be determined a bit later. Given a number j let tj be instant of the corresponding 
Hadamard transformation. We shall obtain the required gate array inserting sequential NOT 
operations on each qubit. Construct a list of Poisson random processes Ap p — 0, 1, . . . , I — 1 
corresponding to qubits and generating time instances for NOT operations on them. But now 
these processes will not be completely independent. At first require that a number of these 
NOT operations on jth qubit preceding tj is even. It will guarantee that a value of jth qubit 
before Hadamard transform will be a'-. Then require that jth and fcth processes generate the 
same time instances in the time frame {tjk^j,k + Atj,k) corresponding to each pair of qubits 
jth and fcth, j > k (qubits are numerated from bottom to top). Call these time frames (j, k)ih 
interval of synchronization. They will have the following properties: 

a) they will not overlap for different pairs, and 

b) each time frame (£?,fc, i/,fc + &tj,k) hes between the instances of fcth and jth Hadamard 
gates, 

c) Atj k ^> 1/ A for a density A of processes A p . 

Appropriately choosing A we can satisfy c). Consider subdivision of intervals between nearest 
Hadamard gates to equal halves and agree that the points of division belong to the left of 
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possible time frames. Agree that (J, k)th interval of synchronization lies in the same half as 
t = (tj + tfc)/2. Then only r < 21 such intervals will lie in the same half and to make them 
not overlapping we only should divide each half to r segments which lengths will be completely 
determined when we choose the values of Atj^, associate each interval with some segment and 
then assume that each of these r intervals belongs to the corresponding segment. Thus we 
ensure conditions a) and b). Now to determine a gate array (look at the picture 3) we need 
only values At,- jfc . 



-e- 



Interval of sinchronization 

-e- 



t 



j,k' 



~e- 



-e- 



h 



tj, k + At 



Picture 3. Gate array implementing QFT . It consists of one-qubit gates 
and different fixed permanent interactions of diagonal form 



At first consider the case when numbers of NOT operators preceding the first NOT in syn- 
chronized interval have the same oddity for jth and kth qubits. Count a phase shift generated 
by this gate array omitting linear summands like Cjd'j or d^b^ which can be easily compensated 
by the appropriate one-qubit gates: Ai,yfc|((l — a^)(l — bk) + a'jbt) + Su n + S neg i. Su n consists 
of linear summands, S neg i issues from interaction outside synchronized intervals and for high 
density A it consists of mutually canceling additions like above. In the case of different oddities 
we obtain a similar expression. In all cases we can choose such values for all At^ (they will 
depend on pairs j, k ) that our gate array fulfils QFT -1 in within linear phase shift and it takes 
a time 0(l 2 ). Composing it with gates for appropriate linear phase shifts before and after this 
gate array we obtain a required implementation of QFT -1 in time 0(l 2 ). 

Given interaction of the form ([1] , A) where p is negative we can reduce this case to the 
considered one by inverting one of interacting gates and adding linear shift. It is straightfor- 
wardly seen that it results in the only change in our constructions as adding linear shift to the 
phase. By the same manner we can implement the direct transformation QFT given a gate 
array implementing its reversal and complexity of the gate array will be about the same. 

Given interactions of the form ([I], B) its application in time frame 1 gives phase shift to 
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pi(l — a'j)(l — bk) +^2(1 — ct'j)bk + P3«j(l — bk) + p^a'jbk that may be reduced by the linear phase 
shift to the above case because p\ + p 4 7^ p 2 + P3- 

At last this technique may be straightforwardly generalized to interactions with Hamiltoni- 
ans that may be diagonalized by one-qubits transformations. Note that this class of interactions 
is not yet the most general at all, say a model of quantum computations with CNOT gates 
cannot be immediately reduced to it. 



5 Simulation physics by means of fixed interaction 



Now we take up an important idea dating back to Feynman ( ||Fej| ): to simulate physics by 
quantum computers. A sketch of such simulating method referring to Coppersmith-Deutsch- 
Shor scheme (see picture 1) for QFT^ 1 was proposed by Zalka ( |Za|1 ) and Wiesner ( |[Wi|| ) . The 
method of QFT implementation presented above gives an easy way for simulation in case of 
linear and quadratic potentials by means of constant and permanent interaction between qubits. 
Hamiltonians with quadratic potentials serve as a good approximation for description of such 
important physical objects as free particle, ensembles of linear harmonic oscillators, free fields, 
complex molecules. Such Hamiltonian for si particles has the form 

where s = 3si is the total number of spatial coordinates qk determining a spatial state of 
system, p^ are impulses and Vjk are constants. We now take up a case k = 1 because the 
general case may be considered similarly. 

At first we remind the main ideas of simulation physics by a quantum computer. We have 
to approximate an action of operator e~ lHt on a wavefunction ipQ where H = H p + H q , H p = 
H q = V{q) ) P = and potential V{q) is a real quadratic function. Without loss of generality 
we can take t — 1. To have a useful approximation we must deal with coordinate or impulse 
basis in the space of state vectors and in both cases Hamiltonian will not be diagonal. In order 
to reduce the problem to the simple diagonal case choose a small time interval At and represent 
our evolutionary operator approximately by 

e~ iH « (e~ iHqAt e - iH P At y/ At . (8) 

Choosing, say coordinate basis we have a diagonal operator H q . Using Fourier transform 
FT : / — ► -7= f^™ e~ vpq f(q) dq and its property to replace derivative d/d q by factor ip we 
can represent an action of impulse part of operator as e~ tHp = FT" 1 e -*p 2A */ 2m px where the 
medium operator has a diagonal form. If we can implement FT and phase shift on —p 2 /2m 
then the sequential application of such operators from (||) gives the required approximation. 

Assume that a wave function ip(q) is defined on a segment (—A, A) and its impulse rep- 
resentation FT t/> is defined on a segment (—B,B). Choosing a small values Aq and Ap we 

2A/Aq 

can approximate it by J2 ip{q a )5 a where S a (q) takes a value 1 on a segment (q a , q a + Aq) and 

a=0 
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zero for other q. Then we can approximate FT by a linear operator whose action on 8 a gives 

2B/Ap 

-k=Aq e tpbqa crb(p) where crJp) is one step impulse function analogous to 5 a . Introducing 

V2n 6=0 

new one step functions for coordinate and impulse by d a (q) = 5 a (q — A), Sb(p) = <Jb(p ~ B) we 
rewrite FT in the form 

d a — -^Aq £ e- i6aA « A % (9) 

that looks similar to QFT. 

Assume that the physical space is grained in coordinate and impulse representations with 
the sizes of grains Aq and Ap. Then our particle may exist only in points of the form q a or 
may have impulse only of the form p^. We associate a position q a ; a = 0, 1, . . . , N = 2 with 
the basic state \a) of I qubits quantum system. For the simplicity choose such units for the 
length that Aq = Ap = ^phxj\[~N and let A = B = ^ttN/2. Then (|) corresponds to QFT of 
the form @ and phase shift to —p 2 At/2m from (H) corresponds to phase shift to — nb 2 At/mN. 
We can implement the both operations by means of fixed interactions because the last has the 
form from above. At last the first operator in @ can be implemented by the same way in the 
case under consideration. 

If we simulate a system with si particles then we should get s% copies of quantum register 
for one particle and fulfill the procedure described above for this joint quantum memory. 



6 Advantages of quantum simulation of wave functions 

A proposed way of quantum solution of Shroedinger equation as well as implementation of 
QFT preserves all advantages of known quantum tricks (look at |f3h| , |Za| , |Wi|| ). Namely, given 
constant interactions between qubits our method of implementation of QFT requires the time 
0(l 2 ). If we restrict ourselves by approximate version of QFT (it may be obtained by omitting 
exponentially small phase shifts) then the correspondent modification of our method takes the 
time 0(1) where a constant depends on the chosen accuracy. 

As for the simulation of wave function the main advantage of quantum method are displayed 
for the case when a simulated system contains many particles. Let we have a system of s one- 
dimensional particles. Its wave function in a fixed time instant has the form ip(xi,X2, ■ ■ ■ ,x s ) 
where Xj denotes a coordinate for jth particle. To store the approximation of this function with 
grain e and arguments bounded by b we need of order iV s bits where N = b/e. At the same time 
quantum method of simulation requires of order log N qubits for each particle and all memory 
will be of order slogiV qubits that is logarithm of classical size. But the proposed method 
shows advantages even in case of one particle. Limit the time frame of the required simulation 
for the simplicity by 1, a coordinate and impulse - by B. To use the formula @ we must have 
only At — ► 0. Then the total time of quantum simulation has an order log 2 (l/e)-^ whereas an 
approximate solution of Shroedinger equation on a classical computer requires scanning of all 
massive of wave function values that is 1/e in each of 1/At passes. Hence quantum methods of 
simulation give almost exponential time and space saving as compared with classical methods. 
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7 Conclusion 



We considered a model of quantum computer controlled only by one-qubit impulses whereas 
interaction between qubits is fixed and remains unchanged in course of computation. Advantage 
of this model is taken of the simplicity of control. It was found a simple way to implement 
Quantum Fourier Transform by such quantum computer. This method makes possible to solve 
Shroedinger equation for linear and quadratic potentials and thus in principle such type of 
computer can simulate systems of harmonic oscillators, free fields and particles and molecules 
consisting of many atoms. 

Further investigations can go to several directions. The first is: to clarify possibilities of 
such simple model of quantum computer, say to find an effective implementation of Grover 
search algorithm ( ||Gr|| ) in the framework of this model (preliminary step in this direction 
was made in | |FG| | but it is not a complete feasible solution). Would be important to extend 
the results to not diagonal Hamiltonians. At last it may be interesting to reformulate some 
areas of quantum mechanics in terms of qubit representation of wave function that we used 
for quantum computing. For example such reformulation assumes the existence of spatial 
grain and traditional difficulty of field theory issuing from the divergence of rows for high 
frequencies would be removed. This reformulation also can turn to be more economical than 
the conventional because the time and space resources required to describe a many particles 
system grows much slower than for conventional case when the number of particles increases. 
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